Negative isotope effect of a BCS-like gap; 
an inelastic light scattering study 
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In a prior study of the boron isotope effect on the transition temperature (T e ) of -RNI2B2C (R = 
Y, Lu) system, a positive isotope effect was observed, e.g., T C ( 10 B) > T c ( n B). BCS theory predicts 
that the superconducting gap at zero temperature (Ao) is proportional to T c . Therefore, the gap 
should also show a positive boron isotope effect. On the contrary, in an inelastic light scattering 
study, we report a negative boron isotope effect on the energy of the gap-like feature in B2 9 
symmetry, e.g., n{ 10 B) < fi( u B). The origin of the effect is discussed in terms of selective probing 
of a non-conventional electron-phonon interaction on part of the Fermi surface, probably near the 
X point in the Brillouin zone. 

PACS numbers: PACS numbers: 78.30.Er; 74.70.Dd; 74.25. Jb 



In a previous study of an isotope effect on the tran- 
sition temperature (T c ) of i?Ni 2 B 2 C (R = Y, Lu) 
system, Q positive isotope exponents ajf = +0.21 ±0.07 
for YNi 2 B 2 C and +0.11 ± 0.05 for LuNi 2 B 2 C were ob- 
served upon changing the mass of the boron site. The 
exponent a^ c is defined as a^ c = — ^"mb ' ^ rela- 
tively large positive value despite the light mass of boron 
is regarded as a good indication that these materials 
are BCS-type superconductors, where phonons involving 
boron atoms play a major role in the pairing mechanism. 

BCS-type theory predicts that the superconducting 
gap in the limit of the absolute zero temperature (A ) 
is proportional to T c . In weak coupling A = 1.76fcsT c , 
where &b is the Boltzmann constant. The gap should also 
show a positive isotope effect upon changing the boron 
mass. We report here the contrary result: In an inelas- 
tic light scattering study, we observe a negative isotope 
effect of the gap-like Raman peak in B 2g symmetry. 

The gap A can be measured by electronic Raman scat- 
tering. The photon scattering process excites two quasi- 
particles out of the condensate at a minimum energy cost 
of 2 A. In the limit of small wave- vector transfer q rela- 
tive to the inverse coherence length and small elastic and 
inelastic scattering rate relative to 2A/h and with ne- 
glect of final state interactions, the Raman spectrum has 
a peak at an energy shift of 2A. Earlier measurements 
on NbSe 2 and A15 compounds showed redistribution of 
the continuum of the Raman spectra and the formation 
of a 2A-like peak upon cooling the samples below T c .^ 
Electronic Raman studies have played important role in 
characterizing the nature of the superconducting gap of 
the high-temper aure superconductors. [0 

Compared with other techniques that measure the gap, 
such as tunneling and photoemission spectroscopy, in- 
elastic visible light scattering has the advantage of mea- 



suring the anisotropy of the gap while being less surface 
sensitive. 

The single crystal samples were grown by the flux- 
growth method j| || and were characterized by resistiv- 
ity, magnetization, and neutron scattering. i?Ni 2 B 2 C 
crystalizes in the tetragonal body-centered space group 
IA/mmm, and phononic Raman analyses have been 
made earlier. Here, we concentrate on the electronic 
Raman scattering from the quasi-particles in normal and 
superconducting states. 

Raman spectra were obtained in a pseudo- 
backscattering geometry using a custom-made sub- 
tractive triple-grating spectrometer designed for very 
small Raman shifts and ultra low intensities. § A 3 
mW beam of 6471 A Kr-ion laser light was focused 
onto a spot of 100 x 100/^m 2 . The temperature rise 
AT of the focal spot above ambient temperature was 
estimated by solving the heat-diffusion equation. AT 
is largest at lowest temperature because of the smaller 
values of the thermal conductivity. || The estimated 
AT is 2.7K at 4K and 0.9K at 14K for YNi 2 B 2 C 
single crystals. The ambient temperature at which the 
Raman continua begin to show the redistribution was 
determined to be 14K, in agreement with this estimate. 
The spectra were corrected for the Bose factor and 
therefore are proportional to the imaginary part of the 
Raman susceptibility. 

The ll B 2g spectra" and "Bi g spectra" presented below 
are actually from the linearly-polarized light scattering 
geometries of XY (B 2g + A 2g ) and X'Y' (B lg + A 2g ), 
respectively. The A 2g contribution to these spectra was 
determined from the sum of XY and XT' spectra minus 
the spectrum obtained in LR (Bi g + B 2g ) circularly- 
polarized scattering geometry. A 2g was found to be neg- 
ligible. 



FIG. 2: Isotope dependance of the the "2A-like" peaks in 
Raman spectra measured at 6 K. a) B^g peaks of B10-Y (thick 
line) and Bll-Y (thin line), b) B29 peaks of BlO-Lu (thick 
line) and Bll-Lu (thin line), c) A\ g peaks of B10-Y (thick 
line) and Bll-Y (thin line), d) B lg peaks of B10-Y (thick 
line) and Bll-Y (thin line). 



no such isotope dependence is observed. On the contrary, 
the Big peak tends to show a positive isotope effect, but 
the greater breadth of the peak hinders determining an 
isotope effect within experimental resolution. 

In our earlier paper, |io[ it was shown that the -B 2g 
peak follows the BCS prediction for the temperature de- 
pendence of the superconducting gap. However, as seen 
in Fig. 2a, it does not follow the BCS prediction for the 
isotope dependence of the gap. Rather, the B^ g peak 
frequency in YNi 2 B 2 C, as well as the leading edge of the 
peak, exhibits a negative isotope effect. There are prior 
reports of some of these experimental results. |l2| 

We found empirically that in the region of the 2 A peaks 
the ratio of Bi g to B\ g Raman intensities can be fit to a 
Lorentzian expression 

l^L = A (1) 

I Bls (uj - u y + T2 [ > 

near the 'gap-like' peak frequency, u>q. 

Figure [| shows the results of fitting the above equation 
to the ratio of the B2 g to B\ g spectra. The fit is excellent 
around the peak and below, deviating much at higher 
frequencies. No such fit is possible for Ib 2 or Ib 1 itself. 
From the fits to several sets of B-2 g /Bi g Raman ratios, 
we obtained the isotope effect exponent due to boron, 
a Raman = _ Q 33 ± q 03 for Y Ni 2 B 2 C and -0.04 ± 0.07 

for LuNi 2 B 2 C, where a § aman = M B is the 

mass of the boron isotopes. Since there is little evidence 
for an isotope effect in the B\ g spectra, we attribute these 
results for a ^ aman primarily to the negative isotope effect 
on the Bi g spectrum. 



FIG. 4: Plot of B2 g peak frequencies (fl) at 6K and the 
onset temperatures of the superconducting transition of six 
different YM2B2C single crystals. The dashed line on the 
plot represents BCS prediction 2Ao = 3.52fc_eT c . 



We have measured the superconducting transition tem- 
peratures using a SQUID magnetometer and eliminated 
the possibility that the particular BIO samples in our 
study have lower T c values than those of Bll samples. 
Figure ^ shows a plot of B 2g peak frequencies (f2) at 
6K and the onset temperatures of the superconducting 
transition of six different YM2B2C single crystals. As 
is clearly seen, Bll-Y samples do have larger B 2g peak 
frequencies than B10-Y samples, in contradiction to the 
BCS prediction. This figure shows that the isotope effect 
of the B 2g peak frequencies is negatively correlated with 
the isotope effect of the T c in YN12B2C. The dashed line 
on the plot represents 2Ao = 3.52ksT c , the BCS predic- 
tion in the weak coupling limit. 

There is no previous report known to us of an observed 
negative isotope effect of the superconducting gap. Neg- 
ative isotope effects have been observed for the transi- 



tion temperatures of the Pd-H(D) system|l3[| and under- 
doped YBa 2 Cu307_5 .Q In the following discussion we 
assume that the B 2g peak is directly due to supercon- 
ducting gap 2A. How, then, can we explain a negative 
isotope effect for the gap measured in B 2g Raman sym- 
metry along with a positive isotope effect for T c and small 
or zero isotope effects for the gap measured in B\ g and 
A\ g symmetries? 

The most likely possibility is that the boron portion 
of the electron-phonon interaction plays out quite differ- 
ently on the various portions of the Fermi surface. It 
would be necessary for the B 2g Raman vertex to probe 
those parts of the Fermi surfaces which have a negative 
boron-isotope effect on the value of the gap. Correspond- 
ingly, the isotope effect would have to be much smaller 
and possibly positive on those parts of the Fermi surface 
probed in A\ g and B\ g Raman geometries. 

There is a simple picture of the electronic structure 
that may give a strong clue to where the B 2g vertex is 
large. According to the band structure calculations for 
LuNi 2 B 2 C]l|, |l(| and YNi 2 B 2 C,|l7| there is a flat band 
that may cross the Fermi surface near the X point [wave- 
vector given by (jr/a, 7r/a, 0)]if there is a small change 
in parameters. This band has strong anti-bonding dda 
interactions as well as d z 2 character on Ni sites. |lq] If 
we take a nearest neighbor, planar, Ni-Ni tight bind- 
ing model with coupling t, we obtain a band disper- 
sion e(k) = 4tcos(fc x a/2) cos(k y a/2). This has a saddle 

point at X. A good starting approximation to the Ra- 
man vertex is the so-called mass approximation, valid 
when the photon energy is less than that of all inter- 
band transitions starting or ending at the Fermi energy. 
It says that the B 2g Raman vertex is proportional to 
d 2 e/dk x dk y . For the assumed band structure, this would 
give d 2 e/dk x dk y = 4i(a/2) 2 sin(A^;a/2) sin(fc y a/2), which 
takes on its maximum value at the X point. 

It would be useful to verify this suggestion with a band 
structure calculation of the mass vertex of the band that 
crosses the Fermi energy near the X point. The next step 
would be to study the coupling between this band and 
the boron phonon modes. There is one such calculation 
in the literature, which emphasizes the role of the Ai g 
zone center boron mode on modulating the angle of the 
B-Ni-B tetrahedra. jOI However, much more needs to 
be done, with a particular look for novel couplings that 
would give a negative isotope effect on the gap associated 
with this band. 

The superconducting gap of borocarbides will not be 
as simple as calculated by BCS theory. The finite inten- 
sity below the gap observed in our Raman measurements 
(Figs. 1 and 2 and pCj]), a nearly T 3 behavior of the 
electronic specific heat,[|l9| and evidence of small gap on 
the small pocket Fermi surface from de Haas- van Alphcn 
measurements |2C| of superconducting i?Ni2B2C (R = Y, 
Lu) system are good indications of the complex nature 
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of the gap of borocarbide superconductors. 

In conclusion, we observed a negative boron isotope ef- 
fect of the gap-like feature in B 2g electronic Raman spec- 
tra from YNi 2 B 2 C single crystals. For LuNi 2 B 2 C, the 
negative isotope effect is weaker, and may be zero within 
our experimental resolution. A likely explanation is that 
the _B 2g Raman vertex selectively probes those portions of 
the Fermi surface which show a negative boron-isotope ef- 
fect of the order parameter in the superconducting state. 
We have suggested that this may occur at the X point, 
the corner of the Brillouin zone in the basal plane. 
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